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Abstract 

The self-interaction-corrected local-spin-density approximation is used to describe the electronic 
structure of dioxides, REO2, and sesquioxides, RE2O3, for the rare earths, RE=Ce, Pr, Nd, Pm, 
Sm, Eu, Gd, Tb, Dy and Ho. The valencies of the rare earth ions are determined from total energy 
minimization. We find Ce, Pr, Tb in their dioxides to have the tetravalent configuration, while for 
all the sesquioxides the trivalent groundstate configuration is found to be the most favourable. The 
calculated lattice constants for these valency configurations are in good agreement with experiment. 
Total energy considerations are exploited to show the link between oxidation and /-electron der- 
ealization, and explain why, among the dioxides, only the Ce02, Pr02, and TbC>2 exist in nature. 
Tetravalent Nd02 is predicted to exist as a metastable phase - unstable towards the formation of 
hexagonal Nd2C>3. 

PACS numbers: 



1 



I. INTRODUCTION 



Even though the rare-earth (RE) elements readily oxidize, they do so with varying 
strength.— Ce metal oxidizes completely to CeC>2 in the presence of air. Pr occurs natu- 
rally as Pr 6 On, exhibiting a slightly oxygen deficient fluorite structure. The stoichiometric 
fluorite structure PrC>2 exists under positive oxygen pressure. The rare-earth oxides from 
Nd onwards, with the exception of Tb, all occur naturally as sesquioxides, RE 2 3 . Tb oxide 
occurs naturally as TbztOy, and transforms into Tb02 under positive oxygen pressure. 

The present study is concerned with the valencies of the oxides of Ce, Pr, Nd, Pm, Sm, 
Eu, Gd, Tb, Dy and Ho. Under suitable conditions, all the RE elements form a sesquioxide, 2 
and there is general agreement that, in the corresponding groundstate, the rare-earth atoms 
are in the trivalent, RE 3+ configuration.^ Each rare-earth atom donates 3 electrons to the 
strongly electronegative O ions, and the remaining 4/ electrons stay strongly localized at the 
rare-earth site. In the lighter lanthanides, the /-electrons are less tightly bound, resulting 
in compounds that display a larger oxygen coordination number, Ce02 and Pr0 2 . Similarly 
in Tb, the extra electron on top of the half filled shell is less tightly bound resulting in a 
valency larger than 3+. 

Only Ce, Pr, and Tb form dioxides, and with respect to the valency of the correspond- 
ing RE-ions the debate is still ongoing. Two conflicting points of view, both based on the 
interpretation of core-level spectroscopy studies, describe the dioxide groundstate as either 
tetravalenl£&£ or intermediate- valent . 8l9ll ° The intermediate- valent interpretation uses the 
spectroscopic data to obtain the model parameters entering the Anderson impurity Hamil- 
tonian.— By solving this many-body Hamiltoniani 2 - the initial state is derived from the 
final-state photoemission spectra. It should be noted however that Wuilloud et a/.— derive 
the tetravalent goundstate for Ce0 2 , using the same method. Marabelli and Wachte r 13 i 14 
conclude from their optical spectroscopy measurements that Ce0 2 is a tetravalent insulator. 
Neutron scattering results on Pr02, by Kern et alM were interpreted in terms of a local- 
ized 4/ 1 groundstate configuration of the Pr ion. Similarly, recent neutron spectroscopy 
measurements of the magnetic exitations in Pr0 2 have also been interpreted in terms of 
a tetravalent 4/ 1 configuration-^ With respect to the XPS results on Pr02, the propo- 
nents of the tetravalent picture explain the data in terms of the coexistence of localized and 
delocalized /-electrons.— 
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The rare-earth oxides find important applications in the catalysis, lighting and electronics 
industries and thus are extremely interesting compounds in their own right. A further reason 
for the strong interest in the RE-dioxides is their relation to the 123-cuprates REB^CusCV 
One notices for example that it is again Ce, Pr and Tb that distinguish themselves from 
the other rare earth elements, as their corresponding cuprates do not support supercon- 
ductivity— Also in this case, no definite picture of the valency of the rare-earth ion has 
yet emerged. Similar to the intermediate-valency interpretation of the PrC>2 data, it has 
been suggested that the suppression of superconductivity in PrBa 2 Cu 3 7 , is associated with 
the existence of intermediate-valent Pr ions.— The determination of the valency configura- 
tion in the dioxides might thus be very useful for the elucidation of the mechanism behind 
superconductivity in the cuprates. 

Any theoretical description of the rare-earth oxides needs to take into account the strong 
on-site /-/ interactions that push the 4 /-electrons towards localization. Using the Anderson 
impurity Hamiltonian, a possible approach is to derive the intraatomic Coulomb energy U 
from experiment, as was done in connection with the interpretation of XPS data.-^ There 
is also a possibility of obtaining the parameters of the Hamiltonian from first principles 
calculations, as was done by McMahan et al.. 19 Electronic structure calculations, based on 
density functional theory have been very successful in describing the cohesive properties of 
solid-state systems with itinerant valence electrons. However the exchange and correlation 
effects of the homogeneous electron gas, which are included in the local-spin-density (LSD) 
approximation, are insufficient to account for the strong correlations experienced by the lo- 
calized 4 /-electrons in the rare-earth compounds. With respect to the rare-earth oxides, the 
band structure calculations that we could find in literature are for CeO-2,-22^ p r 2j 2flJi 
Ce203 , 21l2a , and RE-sesquioxides in general.— From these studies it becomes clear that for 
example Ce02 is best described in terms of itinerant /-electrons, whilst for Ce203, consid- 
ering one Ce /-electron as part of the core leads to better agreement with the experimental 
lattice parameter.— In this case, the /-electron groundstate configuration, and consequently 
the choice of calculational tool, is determined from a comparision to empirical data. Fabris 
et ali^ show that the LDA+U approach gives good agreement with experiment for both Ce 
polymorphs, but the method requires input of a Hubbard U parameter (3 eV for Ce). In 
Pr02 the LSD description reveals a large /-peak at the Fermi level, in disagreement with 
the fact that Pr02 is an insulator. 
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In the self-interaction-corrected (SIC) local-spin-density approac h 26 ! 27 of the present pa- 
per, both localized and delocalized /-electrons are assumed to coexist. The delocalized /- 
levels move in the LSD potential, their exchange and correlation energies are those derived on 
the basis of the homogeneous electron gas, and their ability to form bands leads to a gain in 
hybridization energy. The localized /-levels acquire core-like character by correction of the 
LSD total energy functional for their spurious self-Coulomb and self-exchange-correlation 
energies.— This leads to an additional negative potential term and effectively prevents any 
hybridization. In the SIC-LSD approach, the localized and itinerant /-states are treated 
on an equal footing, and configurations with varying numbers of localized /-states can be 
compared with respect to their total energies. Consequently the preferred groundstate con- 
figuration, as far as the number of localized /-electrons is concerned, can be determined 
from the global total energy minimum. Namely by comparing the total energies resulting 
from different electronic configurations, one can determine whether it is favourable for the 
/-electrons to localize or to contribute to the band formation. The method has previously 
been succesfully applied to describe the valencies of the RE elements^ and compounds (see 
for example [30] and references therein). 

The remainder of this paper is organized as follows. In section II and III respectively the 
results for the RE-dioxides and RE-sesquioxides are presented and compared to experiment. 
In section IV, the oxidation/reduction process is studied, and in section V, we draw our 
conclusions. 

II. THE RE-DIOXIDES 

The SIC-LSD approach has been implemented using the tight-binding linear muffin- 
tin orbital (LMTO) method 31 in the atomic sphere approximation (ASA). The spin-orbit 
interaction is included in the Hamiltonian. Empty spheres are inserted on high symmetry 
interstitial sites. The valence panel includes the 6s, 5p, 5d and 4f orbitals on the rare-earth 
atom, and the 2s and 2p on the O atom, while the higher order i orbitals of O (3d and 
4/), as well as all degrees of freedom of the empty spheres are treated as downfolded.— 
A separate energy panel is used to describe the semicore 5s states of the rare earth. In 
the valence band, we need to distinguish between those /-electrons that are SIC-localized, 
and the delocalized /-electrons, which together with the s, p and d electrons constitute 
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Compound 


E IV — E in 


Vtheo. 




Ce0 2 


-2.40 


39.61 


39.6 


Pr0 2 


-1.44 


39.22 


39.4 


Nd0 2 


-0.65 


39.37 




Pm0 2 


-0.18 


39.15 




Sm0 2 


0.49 


42.19 




Eu0 2 


2.31 


41.87 




Gd0 2 


1.22 


41.08 




Tb0 2 


-0.27 


36.50 


35.6 


Dy0 2 


0.05 


39.65 




Ho0 2 


0.46 


39.06 





TABLE I: Dioxide data: Column 2: Energy difference between tetravalent and trivalent configu- 
rations (in eV). Column 3: Theoretical volume in (A 3 ). Column 4: Experimental volume in (A 3 ) 
from Ref. bj 

the valence electrons. Different electronic configurations will give rise to different nominal 
valencies which in the SIC-LSD approach are defined as: N va i = Z — N core — Ngic, where 
Z is the atomic number, N core is the number of atomic core and semicore electrons, and 
Nsic is the number of localized / electrons. According to this definition, a localized f 1 
configuration of the RE ion will be referred to as trivalent in the case of Ce, tetravalent in 
the case of Pr, pentavalent in the case of Nd, etc.. 

The total energies of all the RE-dioxides from Ce0 2 to Ho0 2 were calculated with a rare 
earth valency assumed to be either 3, 4 or 5. In all these cases the crystal structure was 
taken to be the cubic ffuorite structure. The results are listed in Table I, where the energy 
differences between the tetravalent and trivalent configuration, Eiy — Em, are displayed 
in column 2. Here, a negative value indicates that the compound prefers the tetravalent 
groundstate configuration. The pentavalent configuration is found to be energetically un- 
favourable in all the dioxides. An example of the calculations involved in determining the 
energetically most favourable groundstate configuration is shown in Fig. Q for Pr0 2 . In 
this figure, the total energy of respectively pentavalent (Pr(/ )), tetravalent (Pr(/ 1 )), and 
trivalent (Pr(/ 2 )) Pr0 2 , is plotted as a function of volume. The global energy minimum is 



obtained in the tetravalent configuration. The theoretical equilibrium volume is calculated 
to be Vfeto=39.2 A 3 , which is in excellent agreement with the experimental value V exp =39A 
A 3 . 

A closer examination of the electronic structure reveals why the trivalent configuration 
is energetically unfavourable. In Fig. |21 both the total and /-projected DOS of Pr0 2 are 
shown for the pentavalent (Fig. |2K), tetravalent (Fig. Eb), and trivalent (Fig. Efc) Pr 
configurations. In Fig. with all the /-electrons treated as delocalized, we find the Fermi 
level in the /-peak, in accordance with the LSD calculations by Koelling et a/., 20 and, as 
noted earlier, in disagreement with the observed insulating nature of PrC>2. Localizing one 
/-electron, results in the DOS of Fig. |2t>> where now the Fermi level is situated inside the 
large gap, which forms between the occupied O p states and the unoccupied conduction 
bands of primarily Pr d character. Even though one /-electron has become localized, we 
still observe considerable / hybridization with the O p-band, i. e., the O p-states have large 
spatial extent and their tails reach into the atomic sphere around the RE atom, where they, 
when decomposed into RE-centred spherical harmonics, attain appreciable / character. 

In this respect, both localized and delocalized /-electrons coexist in the tetravalent con- 
figuration, which is in line with the interpretation of X-ray absorption studies™^ Finally, 
in the trivalent scenario of Fig. Et, a further /-electron becomes localized, which results 
in the Fermi level moving down into the p-band. The O p-band states are depopulated to 
facilitate the formation of Pr f 2 ions, i.e. charge transfer is imposed on the system, and the 
associated cost in energy is significantly larger than the gain in /-localization energy. The 
calculated energy minimum for the trivalent f 2 configuration is situated more than 1 eV 
above the tetravalent energy minimum (E/y-E///=-1.44 eV), implying that localizing one 
more /-electron with respect to the tetravalent configuration results in a loss of hybridiza- 
tion and electrostatic energy (~ 2.5 eV) that considerably exceeds the corresponding gain 
in SIC localization energy (~ 0.95 eV per /-electron). According to the proponents of the 
intermediate valency, the groundstate of Pr02 consists of a 50:50 mixture of 4/ 1 and 4/ 2 
states , with an average 4/ occupancy of approximately 1.5 electrons. In the present 
picture this would imply the trivalent and tetravalent configurations to be nearly energeti- 
cally degenerate. This does not occur, as our total energy results point towards a relatively 
solid tetravalent groundstate configuration for Pr02- The calculated band gap for Pr02 
is approximately 1.1 eV, considerably larger than the 0.262 eV derived from conductivity 
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FIG. 1: Total energies (in eV per formula unit) as a function of volume (in A 3 ) for PrC>2, assuming 

the Pr ions to be pentavalent (/ , dotted line), tetravalent (f 1 , solid line), and trivalent (f 2 , dashed 

line), respectively. The global energy minimum is obtained in the tetravalent configuration. 
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FIG. 2: Total DOS (solid line) and /-projected DOS (dotted line) for Pr02, with the Pr ions in a) 
the pentavalent (/ ) configuration b) the tetravalent (f 1 ) configuration and c) the trivalent (/ 2 ) 
configuration. The energy is in units of eV, with zero marking the Fermi level. 

measurements by Gardiner et aL— 

For Ce02 we find a clearly preferred tetravalent groundstate configuration (the trend 
towards derealization of the /-electrons is even stronger than in PrO^), i.e. Ce02 is best 
described in the LSD approximation, in line with results from earlier band structure cal- 
culations.— >2i The calculated volume in this configuration is in good agreement with the 
experimental value, as can be seen from Table I. 

In conclusion, our calculations do not confirm an intermediate-valent groundstate for 
either PrC>2 and CeO^. Koelling et alM argued that, the intermediate valency scenario, is 
related to an ionic description of the RE-oxides, which can not account for the covalent f-p 
bonding. In the ionic picture, valency is denned as the number of electrons that have been 
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transferred from the RE-atoms to the O-atoms, i.e. /-electrons do not participate in the 
bonding, and can only exist as localized at the RE-sites. A given configuration thus has an 
integral number of /-electrons, and a non-integral number of /-electrons can consequently 
only result from an intermediate-valency scenario. In the SIC picture both localized and 
delocalized /-states coexist. The delocalized /-electrons are allowed to participate in the 
band formation, and occur as part of the tails of O p states. The overall number of /-electrons 
is non-integral, in analogy to the intermediate valent ionic picture, but as a consequence of 
the p — f mixing. 

Apart from Ce and Pr, Tb also forms a dioxide, albeit under positive oxygen pressure. 
We find the tetravalent groundstate configuration also to be energetically most favourable 
for Tb02, although here the energy difference, between the trivalent and tetravalent config- 
urations is less pronounced than for Pr02- Again the volume, calculated for the tetravalent 
configuration is in good agreement with the experimental value. The remaining RE dioxides 
do not form in nature, and the determination of their groundstate configuration might seem 
of academic interest only. The energy differences, Ejy-E/jj, in Table I, indicate that Nd02, 
and PH1O2 would also prefer the RE 4+ ion configuration. One notices however a decreasing 
affinity for the tetravalent configuration from Ce02 to Pm02, and from Sm02 the energy 
difference changes sign, indicating that now the trivalent configuration is energetically more 
favourable. With increasing nuclear charge, the /-electrons become more tightly bound to 
the RE atom, and the decreased overlap with neigbouring atoms results in a reduced gain 
in binding energy. Eventually it becomes more favourable to localize an extra electron, 
and gain the corresponding SIC energy, which is what happens in Sm02, with a resulting 
trivalent groundstate configuration. In the late RE dioxides, only Tb0 2 is observed to be 
tetravalent, which is due to the fact that the Tb ion prefers a half- filled /-shell. 

III. THE RE-SESQUIOXIDES 

Below 2000 C, the rare-earth sesquioxides adopt three different structure types.— The light 
RE crystallize in the hexagonal La203 structure (A-type), and the heavy RE crystallize in 
the cubic Mn203 structure (C-type), also known as the bixbyite structure. The middle RE 
can be found in either the C-type structure, or the B-type, which is a monoclinic distortion 
of the C-type structure. Conversions between the different structure types are induced under 
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Compound 


Eiv ~ E in (eV) 


Vhexag. 


(A3) 






hexagonal 


Theory Expt. 


Theory Expt. 


Ce 2 3 


0.38 


76.4 


79.4 


87.27 


87.0 


Pr 2 3 


0.78 


75.6 


77.5 


85.08 


86.7 


Nd 2 3 


0.94 


74.0 


76.0 


83.52 


85.0 


Pm 2 3 


0.97 


72.9 


74.5 


82.92 


83.0 


Sm 2 3 


1.09 


72.0 




81.27 


81.7 


Eu 2 3 


1.13 


70.5 






80.2 


Gd 2 3 


1.29 


68.8 




79.43 


79.0 


Tb 2 3 


1.12 


67.6 




78.21 


77.2 


Dy 2 3 


1.21 


66.3 




77.78 


75.9 


Ho 2 3 


1.36 


65.2 




79.72 


74.7 



TABLE II: Sesquioxide data: Column 2: Energy difference between tetravalent and trivalent con- 
figurations (A-type structure). The experimental c/a ratio was used in the calculations. Columns 
3 and 4: A-type theoretical and experimental volumes. Columns 5 and 6: C-type theoretical and 
experimental volumes. Experimental volumes are from j^ . 

specific temperature and pressure conditions.— We investigated the electronic structure of 
both the A-type and C-type sesquioxides. The C-type bixbyite structure was approximated 
by the ffuorite RE0 2 structure, with 1/4 of the O atoms removed. 

Starting with the hexagonal A-type structure, we see from Table II, column 2, that all the 
sesquioxides prefer the trivalent groundstate configuration, in agreement with experiment. 
The degree of trivalency, E/y-E///, increases from Ce 2 3 to Gd 2 3 , then decreases slightly 
at Tb 2 3 , to inrease again through Dy 2 3 to Ho 2 3 . Similar trends were obtained for 
the cubic C-type sesquioxides. The fact that the energy difference is less for Tb 2 3 than 
Gd 2 3 , indicates that the tetravalent half-filled /-shell configuration becomes relatively more 
important. For this same reason, trivalent Gd 2 3 is energetically very favourable. 

Apart from the extraordinary stability of the half-filled shell configuration, the general 
tendency towards trivalency is clearly related to the increasing localization of the /-electrons 
with atomic number. The highly directional /-orbitals are only partially able to screen each 
other from the attractive force of the nucleus, which results in a steadily increasing effective 



in 



nuclear charge with increasing number of /-electrons. The increase in localization leads to 
the well known lanthanide contraction, i.e. the decrease in ionic radii across the rare earth 
series, which is also reflected in the lattice parameters of the sesquioxides. In Table II, 
the experimental and theoretical volumes are compared, for the A-type (columns 3 and 4) 
and C-type (columns 5 and 6) sesquioxides. The data are also illustrated in Fig. 3. Good 
agreement is seen with respect to both the overall trends and the absolute values. For the 
A-type hexagonal structures, the experimental c/a ratio was used, i.e., no optimization of 
this parameter was attempted. From Fig. |3J we notice that the agreement between theory 
and experiment is considerably better for the C-type structure than for the A-type one. 
This might be related to the ASA used in the calculations. This geometrical approximation 
is likely less reliable when applied to the hexagonal A-structures, than when applied to the 
higher symmetry cubic C-structure. Figure |3] also includes the equilibrium volumes calcu- 
lateded by Hirosaki et a/.—, using the projector augmented-wave (PAW) pseudopotential 
method, with the localized partly filled /-shell being treated as part of the core. For the 
A-type structure, the SIC-LSD calculations consistently yield lower volumes than observed, 
while the PAW values are consistently above. For the C-type structure, the present SIC-LSD 
calculations agree somewhat better with observations than the PAW results, in particular 
for the earlier REs, which may be due to the core approximation made for the /-electrons 
in the work of Ref. [25]. and which may be too restrictive in the earlier RE sesquioxides. 

The DOS and bandstructures for all the sesquioxides are quite similar and as a repre- 
sentative example in the Figs. |U and El we present those of Nd203, as calculated in the 
trivalent ground state 4/ 3 configuration. The broad band below the Fermi level originates 
from the O-p states. Hybridization and charge transfer result in this band being completely 
filled, leaving behind an empty /-peak situated in the gap between the valence and (non- 
/)conduction bands. We find all the sesquioxides to be insulators. 

In Table we compare the calculated band gaps for the A-type structure, to the ex- 
perimental values, as obtained respectively from optical (column 4) and conductivity^^ 
(columns 5 and 6) measurements. The two sets of conductivity measurements refer re- 
spectively to temperature regions above 3 ^ (column 5) and below 37 (column 6) T=800 K. 
We notice a considerable discrepancy between the optical and conductivity measurements, 
and the two sets of results have been interpreted slightly differently. In Fig. El we show 
schematically the various transitions that can be envisaged to occur. The actual transition 
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FIG. 3: Calculated equilibrium volumes of the rare-earth sesquioxides, crystallizing in the hexag- 
onal A-type structure (dotted line), and the cubic C-type structure (full line). The stars, circles, 
and squares refer to the PAW, 25 SIC-LSD (present), and experimental results^ 

depends crucially on the position of the occupied and empty /-states, respectively /" and 
/ n+1 . The interpretation of the optical data is that the empty /-levels are situated above 
the conduction band minimum (CBM), and the energy gap E g is entirely determined by the 
gap between the occupied /- and the CBM. If the /-states move below the valence band 
maximum (VBM), the band gap becomes largest (i.e. v — > c). The activation energies 
E a in column 5, obtained from conductivity measurements (using conductivity~ e -E a /2kT^ 
indicate that, for some compounds, the transitions is from (v — > /), i.e. that the empty 
levels can be situated in the gap beteen VBM and CBM. 

One could speculate whether the absence of (v —>■ /) transitions in the optical spectrum 
is due to this transition being optically forbidden, whilst no such constraint exists for the 
thermal excitations, which go into the unoccupied /-bands. However this can not explain 
the relatively large difference in the gap energy, observed for example in the case of Gd203, 
where both optical and conductivity mesurements agree on the nature of the gap being 
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FIG. 4: Total DOS of NCI2O3 in the trivalent (/ 3 ) configuration in the A-type hexagonal structure 
at the calculated equilibrium volume. Energy is in eV, with marking the mid-gap position, while 
DOS is in states per formula units and eV. 




FIG. 5: Bandstructure of A-type Nd203 in the trivalent configuration. 

(v —>■ c). From the theoretical data we see that the gap between valence and conduction 
bands (v — > c) shows a rather slow increase across the RE series, whereas the gap between 
valence and /-band {v — > /) decreases from Ce 2 03 to EU2O3, and from Gd2C>3 to DV2O3. For 
Ce203, Pr 2 03 and Gd2C>3 no separate edge for the non-f conduction bands was found in the 
calculations. The direct comparison between calculated and experimental gaps is tricky. The 
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FIG. 6: Schematic representation of the RE sesquioxide band gap. The f n and f n depict 
respectively the occupied and empty /-bands. 



SIC-LSD approach cannot accurately determine the removal energies of localized states, and 
the bare /-bands always appear at too high binding energies due to the neglect of screening 
and relaxation effects.— With the SI corrected /-states situated below the valence band 
edge, the calculated band-gaps are of either v-f or v-c character. The observed trend can be 
explained by the fact that with increasing atomic number the /-orbitals become increasingly 
localized, and the energy of the /-shell decreases. At Gd2C>3 all seven majority spins are 
Si-corrected, and the minority states are situated above the conduction band minimum. The 
v — ► / gap starts decreasing again from Tb 2 3 as the miniority states are pulled into the 
gap between valence and conduction band. 

Our calculated values for v c are consistently smaller than the optical E g . If the 
optical E s 's turn out to be a better representation of the actual band gaps, the reason for the 
difference with our calculated results could be twofold: The LSD is known to underestimate 
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Compound 








E gap (eV) 






Theory 




Experiment 






«-/ 


v — > c 


optical 


a (T=800-1200 K) 


cr (T< 800 K) 


Ce 2 3 


3.20 


3.20 


2.4 






Pr 2 3 


2.62 


2.62 


4.6 




0.8 C 


Nd 2 3 


2.13 


3.58 


4.4 


2.36 A (/ c) 




Pm 2 3 


1.66 


3.77 








Sm 2 3 


1.04 


3.80 


5.0 


2.12 B (v -> /) 


1.2° 


Eu 2 3 


0.63 


3.95 


4.4 


1.84° (v -► /) 


1.2° 


Gd 2 3 


3.13 


3.13 


5.3 


2.64° (v -► c) 


1.0° 


Tb 2 3 


2.27 


4.32 


4.1 


2.24 s (/ - c) 


0.8° 


Dy 2 3 


1.71 


4.39 


4.9 


2.82 c (« -► c) 




Ho 2 3 


1.26 


4.48 


5.2 


2.74 c (v -► c) 


1A C 



TABLE III: Sesquioxide band-gaps. Theoretical values for A-type structure: valence to /-band 
(column 2), valence to conduction band (column 3). Experimental values: optical measurements 
(column 4), columns 5 and 6 from conductivity measurements above and below 800 K respectively. 

gaps between occupied and unoccupied bands, and in addition, the unoccupied / bands may 
be positioned extraordinarily low due to the neglect of correlation effects beyond the LSD 
for these. In order to discuss in more detail the accuracy of the calculated energy gaps, one 
would need to further investigate the reason behind the considerable discrepancy between 
the optical and conductivity measurements. 

IV. OXIDATION 

So far, the present theory has demonstrated that all the sesquioxides prefer the trivalent 
groundstate configuration, whilst the dioxides can be separated into tetravalent light RE- 
dioxides, and trivalent heavy RE-dioxides, with the exception of Tb0 2 , which also prefers 
the tetravalent groundstate. From experiment however we know that all the RE-sesquioxides 
are found in nature, whilst the only dioxides that occur naturally are Ce0 2 , Pr0 2 and Tb0 2 . 
With respect to the valencies, these two sets of information lead us to conclude that the 
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RE-sesquioxides are trivalent, and the naturally occurring RE-dioxides are tetravalent, i. e., 
the oxidation process from sesquioxide to dioxide goes hand in hand with the derealization 
of an extra /-electron. The extraordinary catalytic properties of Ceria are according to 
Skorodumova et alM based on the coupling between /-electron localization and oxidation. 
Recent calculations^ have similarly shown that, the possible existence of the higher oxidation 
state, PUO2+X,— can be explained by the occurence of a Pu pentavalent state, which is 
reduced to a tetravalent state in PuC>2. The question remains however, why the dioxide 
only exists for Ce, Pr, and Tb. For example, from the data in Table I it follows that Nd in 
NdC>2 also readily takes the tetravalent configuration, but nevertheless this compound does 
not occur naturally. 

To discuss this issue, we consider the oxidation process: 

RE 2 3 + |0 2 # 2RE0 2 . (1) 

The balance of this reaction in general will depend on the Gibb's free energy of the reactants 
at the given temperature and pressure. The ab-initio calculations of these quantities are 
beyond the capability of the present theory. However we can still to some extent analyze 
the reaction by looking at the zero temperature and zero pressure limit. In that case 
the free energy difference between the reactants reduces to the corresponding total energy 
difference, as obtained by the SIC calculations: 

E ox = 2E SIC (RE0 2 ) - E SIC (RE 2 3 ) - // . (2) 

Still, the comparison of the dioxide and sesquioxide total energies is not straightforward due 
to the approximations of the present SIC-LSD implementation, most notably the ASA. The 
total energies depend on the ASA radii chosen for the constituent atoms. The problem is 
most severe for the low symmetry hexagonal structure of the sesquioxides, and is the most 
probable reason why the calculated equilibrium volumes of the hexagonal sesquioxides are 
in somewhat poorer agreement with experiment, cf. Table II. Unfortunately, the SIC-LSD 
approach has only been implemented within the tight-binding LMTO-ASA method and not 
in its full potential version FP-LMTO, where the problems associated with the ASA would 
not occur. However, in the /° groundstate configuration, the relative error introduced by the 
ASA, can be estimated from a comparision of the total energies, as obtained respectively with 
a FP-LMTO method and our LMTO-ASA implementation of the SIC-LSD method. This 
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REO2/RE2O3 


E o:r (A-RE 2 3 ) 


E OX (C-RE 2 3 ) 


Ce0 2 /Ce 2 03 


-1.90 


-3.54 


Pr0 2 /Pr 2 3 


-0.14 


-1.90 


TV T 1 S~\ /TVT 1 1 — V 

Nd0 2 /Nd 2 3 


0.54 


-0.54 


Pm0 2 /Pm 2 3 


0.00 


0.27 


DmU 2 /bm 2 (J3 


-0.82 


0.68 


Gd0 2 /Gd 2 3 


-10.88 


0.14 


Tb0 2 /Tb 2 3 


-8.16 


-0.27 


Dy0 2 /Dy 2 3 


-4.08 


0.27 


Ho0 2 /Ho 2 3 


-0.68 


0.00 



TABLE IV: Oxidation energies in eV (from equation (2)). Column 2: assuming the oxidation 
process involves the A-type sesquioxide, Column 3: assuming the oxidation process involves the C- 
type sesquioxide. Negative/positive energies indicate that the dioxide/sesquioxide is energetically 
most favourable. The bold characters indicate which type of sesquioxide will preferentially be 
formed, if at all. 

provides a FP-ASA total energy correction, which can be used to calibrate the corresponding 
energies calculated with the SIC-LSD. Since the energy differences between the valency 
configurations of a given compound are almost unaffected by the chosen ASA radii, for a 
given compound the same correction can be applied to the total energies in the different 
localization configurations. 

In the following, we assume that the dioxides and the C-type sesquioxides, which both 
crystallize in a cubic structure, are adquately described within the ASA approximation. The 
FP-ASA correction is therefore only applied to the A-type sesquioxide, which crystallizes in 
the hexagonal structure. In Table ITVl the resulting oxidation energies (J2J) assuming either the 
A or C-type sesquioxides are given in columns 2 and 3 respectively. The chemical potential 
used, is no = —6.12 eV (relative to free atoms). A negative /positive energy balance means 
that the formation of the dioxide/sesquioxide is energetically most favourable. From the 
table we see that in the case of Ce, Pr, and Tb, the dioxide is preferred energetically with 
respect to both the A- and C-type sesquioxide. We also find that for Ce and Pr the A-type 
sesquioxide is least unfavourable with respect to the dioxide (as indicated by the relatively 



17 



smaller oxidation energy), whilst for Tb it is the C-type sesquioxide that is closest in energy 
to the dioxide. This is relevant since it indicates which sesquioxide structure will form in 
the reduction process. For Nd, we find the dioxide to be more favourable than the C-type 
sesquoxide, but less favourable than the A-type sesquioxide, i.e. at T=0, Nd exists as A-type 
Nd203. Nd02 may exist as a metastable phase, since steric effects may hinder the reduction 
of this compound into the hexagonal sesquioxide. As discussed, up to now the synthesis of 
Nd02 has not been successful. For the remaining RE, the C-type sesquioxide is more stable 
than the dioxide, as indicated by the positive oxidation energies. The relevant oxidation 
energies (indicated by boldface in Table HVjl are plotted in Fig. as a function of rare earth 
ion. The overall picture is that Ce, Pr, and Tb exist as dioxides, whilst the remaining RE 
prefer the sesquioxide, in agreement with experiment. It also shows that with respect to 
the sesquioxides, the later RE from Pm onwards crystallize in the C-type structure, whilst 
reduction of the early REO2 will result in sesquioxide crystallizing in the A-type, which 
again would be in good agreement with experiment. 

There is some uncertainty in our results associated with the value of the oxygen chemical 
potential. We have conveniently used a value of jio = —6.12 eV, as it allows us to indicate 
the trend in oxidation energies more clearly. Actually, using a FP-LMTO calculation, with 
the O2 molecule put into a lattice with large lattice constant, we find \io = —4.76 eV, i.e. 
1.36 eV less than the value used in Table ITVl Using this value in Fig. Ogives a qualitatively 
slightly different picture, as the oxide/sesquioxide borderline is now situated 1.36 eV higher 
(solid line), which means that in this scenario, the dioxide is always preferred. Again other 
LSD total energy calculations give a binding energy of 7.48 eV for the O2 molecule, and 
the experimental value is 5.17 eV.— . We should also notice here, that the calculated total 
energies apply to T = 0, and under actual conditions, the reduction (left side of equation (1)) 
will be favoured by the higher entropy, as it leads to release of O2. Higher temperature thus 
will work in favour of sesquioxide formation. It seems that the existence of tetravalent Ce02, 
Pr02 and Tb02 is an indication that their sesquioxides will readily oxidize. Conversely, we 
could interpret the trivalent configuration of the other RE-dioxides as evidence that their 
sesquioxides will not further oxidize. 
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V. CONCLUSION 



In conclusion, we used the SIC-LSD method to analyze the valencies of the rare-earth 
dioxides and sesquioxides. We find from total energy considerations, that CeC>2, PrC>2, and 
Tb0 2 prefer the tetravalent configuration, with respectively zero, one, and seven localized 
/-electrons. The intermediate valency scenario is not found to be energetically viable, but 
non-integral total /-occupancies do occur due to covalent p — f mixing. The sesquioxides 
are all found to prefer the trivalent configuration, leading to insulating groundstates in 
agreement with experiment. The calculated equilibrium volumes for both the dioxides and 
the sesquioxides are in good agreement with the experimental values. The oxidation reaction 
from sesquioxide to dioxide is found to be accompanied by /-electron derealization, and 
consequently only occurs for the early lanthanides and Tb, for which the /-electrons are 
least tightly bound. 
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FIG. 7: Oxidation energies, E ox in Eq. iffl). for the rare earths Ce to Ho. The full circle is for the 
cubic C-type sesquioxide, the empty circle for the A-type hexagonal structure of the sesquioxide. 
Negative values indicate that the dioxide is stable. The chemical potential of free O is taken as 
fio = —6.12 eV. The solid line just below 1.5 eV indicates the dioxide/sesquioxide borderline in 
case the FP-LMTO calulated value [io = —4.76 eV is used. 
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